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order  to  assess  the  possible  utility  of  gliding  parachutes,  It  was 
necessary  to  develop  a  set  of  methods  for  predicting  behavior  of  gliding 
parachutes..  Gliding  parachutes  presented  some  unique  problems  not  covered 
by  the  technology  base  existing  In  1968.  At  that  time  a  few  brave  jumpers 
were  experimenting  with  squares.  The  Army  was  encountering  difficulty  in 
obtaining  consistent  performance  from  a  small  gliding  system  for  cargo.  On 
some  days  it  worked;  on  others'  it  did  not.  Engineers  with  parachute  exper¬ 
ience  viewed  the  problem  mainly  as  inconsistent  deployment.  Engineers  with 
aircraft  experience  thought  the  problem  was  poor  navigation  and  lack  of 
power.  Neither  field  of  engineering  was  wholly  appropriate. 


Based  on  the  observations  that  canopy  opening  was  not  too  different 
from  conventional  parachutes  in  most  critical  respects  and  that  canopies 
tended  to  glide  with  reasonable  steadiness,  we  chose  to  concentrate  first 
on  the  problem  of  guidance.  As  a  first  step,  we  worked  out  a  graphical 
(manual)  method  of  plotting  trajectory  points  for  a  gliding  canopy  drifting 
with  the  wind  while  homing  on  a  target.  A  visiting  student  helped  with  the 
manual  labor.  Soon  we  transferred  the  technique  to  a  programmed  calculator 
though  graphs  still  had  to  be  drawn  by  hand.  The  technique  was  crude,  but 
it  effectively  demonstrated  some  key  features  of  radial  homing  as  applied 
to  low-speed  gliding  parachutes .(1)  This  was  our  first  use  of  "computer" 
simulation.  Later  we  moved  to  a  fancier  programmable  calculator  with  an 
X-Y  plotter  and  studied  more  exotic  guidance  schemes .(2) (3) 


However,  as  we  investigated  guidance  schemes  it  became  obvious  that  we 
needed  more  knowledge  of  the  flight  dynamics  of  gliding  parachutes.  Pub- 
lished  test  data  was  inadequate.  A  test  program  appeared  to  be  too  expen- 
O—  sive  and  time-consuming.  Also,  we  were  not  sure  what  parameters  should  be 

CP  measured  in  testing  -  nor  how  to  measure  them.  Thus,  we  began  the  develop- 

C-P  ment  of  analytical  and  computational  techniques  for  study  of  all  aspects  of 

.  gliding  canopies.  This  development  effort  has  included  models  for  study  of 

L  j  longitudinal  stability  and  for  study  of  unrestricted  flight  in  six-degrees- 
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*  jO.f-f i reedom  (6D0P)  which  have  been  verif led  recently  by  experiments!  i 

Stability  Analysis 

\  The  longitudinal  stability  is  the  key  aspect  determining  glide  per¬ 
formance  of  gliding  parachute  systems. (5)  A  balance  most  be  achieved  be¬ 
tween  the  various  momenta  siTthat  the  parachute  maintains  a  steady  glide  at 
‘an  angle  of  attack  permitting  useful  glide  velocities  aind  returns  to  this 
angle  if  disturbed.  The  factors  that  determine  the  condition  of  balance 
or  stability  for  a  parachute  are:  i 

1.  The  orientation  and'  ldcatioh~of~lift  and  drag  forces  acting  on  the 

canopy;  ! 

2.  The  magnitude  of  materiaf  mass_ahd  its  distribution  (Including 
payload  mass) ; 

3.  The  magnitude  of  entrained  mass  within  the  canopy  which  is  offset 

by  buoyancy  in  steady  flight  but  nust  be  include*i~in ~deteralnlng~tKe~ Idea-: 
'tion  of  the  system  center  of  mass;  { 

A.  The  location,  magnitude  and  orientation  of  payload  drag.  I 

(The  parameters  available  which  can  be  adjusted  to  achieve  stability  or 
"trim"  at  a  desired  glide  condition  differ  significantly  from  those  avail-  : 
able  in  aircraft  design.  The  key  parameter  is  the  aerodynamic  force  on  the 
^canopy.  For  the  Parafoil  canopy  (Fig  1),  wind  tunnel  data  from  nearly  full 
jscale  canopies  has  been  publlshed(4)  showing  the  variation  of  lift  and  drag 
components  of  aerodynamic  force  with  angle  of  attack.  A  pair  of  empirical  ! 
curves  was  fitted  to  the  data  to  allow  smooth  calculation  of  forces  at  any  I 
angle  of  attack  (Fig  2).  Combined  with  an  analysis  of  geometry  and  mass  j 
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d istribution,  this  permits  determination 
of  stable  conditions  with  variation  of 
control  parameters. 

There  are  two  control  parameters 
which  may  be  varied  independently  or  in 
combination  J:q  achieye_stabillty.  One  I 
is  the  angle  of  incidence  (  T )  controlled 
by  basic  suspension  geometry  holding  the 
chord  plane  of  the  canopy  at  an  angle  to  j 
the  axis  connecting  the  canopy  and  pay- 
load  centers  of  mass.  It  Is  convenient 
to  define  this  angle  as  the  inclination 
of  the  chordllne  below  a  line  perpendicu¬ 
lar  to  the  axis  passing  between  the _ 

centers  of  mass.  For  a  particular  i 

setting  of  incidence  angle,  longitudinal  ! 

r.  -  -  static  stability  is  indicated  by  passage  ! 

and  P tag  Cot^icCwtA  of  the  moment  coefficient  (Cy)  through  j 

aero  with  a  negative  slope  against  angle  j 
of  attack  (Fig  3) .  The  second  control 
parameter  is  the  flap  or  trailing  edge  I 
deflection  which  alters  the  lift  and 
drag  forces  produced  at  any  given  angle  | 
of  attack.  Comparison  of  the  moment  co-  j 
efficient  variation  in  the  two  cases  of  I 
zero  and  full  flap  (Fig  4)  shows  how 
variations  of  either  or  both  of  these 
parameters  can  affect  the  value  of  the  j 
stable  angle  of  attack.  Sport  canopies  i 
are  normally  designed  with  set  angles  of 
incidence  and  with  flap  deflection 
varied  in  flight  to  adjust  the  steepness 
of  the  glide.  We  have  tested  canopies 
designed  with  variable  suspension  geo¬ 
metries  for  in-flight  variations  of  angle 

.  of  incidence. 

Static  Ztahltity 

(Wo  F te  jJ  The  steady-state  glide  velocities 

attainable  at  the  various  stable  angles 
of  attack  (following  some  oscillation)  show  that  a  wide  range  of  perform¬ 
ance  is  available  using  variable  incidence  (Fig  5  and  6).  Optimum  per¬ 
formance  for  a  given  canopy/payload  combination  can  be  obtained  by  adjust¬ 
ment  of  both  control  parameters. 


Distance  separating  the  payload  and  canopy  mainly  affects  the  nteep- 
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Fig  4  Static  Stabitlty 
[Full.  Flap) 


ness  of  the  . slope  .of  _the .pitching .  monent _ 

with  little  effect  on  the  trim  angle  of 
attack.  Payload  drag  haa  a  moderately  sig¬ 
nificant  effect.  In  this  case  a  partially 
streamlined  load  Is  assumed. 

- The  stability  .analysis  has  been  used 

throughout  the  flight  slanlation  studies 
Including  the  6D0F  studies  discussed  below 
and  has  been  verified  In  the  experimental 

-studies -below. _ In  addition  to  Indicating 

the-paraaeter  settings  needed  to  attain  any 
desired  glide  performance,  the  stability 
analysis  has  Indicated  two  potential  problem 
areas  requiring  further  study  in  dynamic 
s loulatlons .and  flight  tests. 

1.  Though  larger  systems  have  the  same; 
wing  loading  as  small  systems  In  order  to 
retain  similar  glide  velocities,  canopy  air 
mass  varies  as  the  cube  of  the  linear  scale.! 
Thus  large  systems  will  have  different 
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static  and  dynamic  stability  which  must  be 
considered  in  their  design. 

2.  Unlike  aircraft  which  start  their 
flight  under  well-controlled  conditions, 
gliding  parachutes  start  flying  at  condi¬ 
tions  differing  greatly  from  their  design 
trim  conditions.  Bence,  undesirable  bi¬ 
stable  conditions  may  be  encountered 
Inhibiting  attainment  of  the  desired  steady 
glide  condition  unless  parameters  are  care¬ 
fully  set  or  some  means  of  adjustment  is 
provided . 

6 DOF  Simulation 


Fig  5 
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A  six-degree-of-freedom  (6D0F)  flight 
simulation  program  has  evolved  out  of 
earlier,  more  limited  simulation  studies  and 
theoretical  analyses.  The  key  features  and 
theoretical  bases  of  the  program  are  dis¬ 
cussed  below.  The  studies  using  the  pro- 
have  yielded  new  Information  as  to  the  nature  of  dynamic  stability, 
control  effectiveness,  and  effectiveness  of  guidance  schemes .  The 
six  degrees  of  freedom  are  the  three  translations  and  the  three  rotations 
of  the  system  of  two  "rigidly  connected"  bodies . 
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Fig  6  Range  o£  Equilibrium 
Velocity  [Full  Flap) 
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Computer  Program 

The  6DOP  simulation  program  for 
gliding  parachutes  features  scaling 
functions  for  geometry  and  mass  dis-  i 
tribution,  allow8_selection  of  pay- 
loads  and  associated  drag  areas, 
selection  of  variable  and  gustlng  two- 
dimensional  wind  functions,  presents  a 
.deployment-window  plot  for  selection  ; 
..of  Initial  conditions,  features  sev¬ 
eral  pictorial  display  options  during  ; 
flight  computation,  includes  a  variety 
of  manual  and  automatic-control  .op- . 
_tiona,  and  presents  plots  of  the  para¬ 
meter  histories  in  a  flexible  format.  ] 
A  recent  modification  allows  simula¬ 
tion  of  random  errors  on  measured  | 
parameters  used  in  guidance  computa-  i 
tlons . 

I 

i 

Applications  of  the  program  in-  I 
elude  study  of  flight  mechanics  and 
study  of  guidance  techniques.  For  j 
example,  to  study  the  turn  control  re-j 
sponse  of  a  particular  canopy/payload  ' 
system,  calm  air  is  selected,  manual  ! 
control  mode  Is  selected  and  a  pictor¬ 
ial  viewpoint  fixed  In  space  near  the  ; 
expected  flight  path  Is  established. 
Generally,  a  view  is  selected  showing 
a  line  drawing  of  canopy  and  payload  ' 
relative  to  an  inertial  background 
(Fig  7) .  The  display  can  be  either  a  ; 
centered  single  image  against  a  moving 
background  or  of  multiple  frozen  lm-  j 
ages  passing  across  a  fixed  back-  i 
ground .  The  preselected  type  of  con-  | 
j trol  deflection  is  Input  via  cross  hairs  against  a  calibration  scale,  with 
i Inputs  possible  at  preselected  Intervals  during  the  flight.  The  rate  of  | 

>  change  of  the  control  deflection  Is  controlled  by  a  time-dependent  func-  j 
tion.  The  operator  observes  the  progress,  sets  desired  control  deflectlonq 
and  changes  the  view  point  or  format  as  needed.  Hard  copies  of  the  dis¬ 
play  are  made  as  required.  Values  of  key  data  parameters  are  also  dis¬ 
played  with  the  pictorial  view  in  one  of  several  formats.  When  the 
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naneuver  jlB .  completed  the  operator _ ! 

a electa  the  output  mode  in  which  a 
ground  track  and  hl8torlea  of  any  of  , 
■  31  key  parameters  are  plotted.  The  , 

operator  selects  certain  combinations 
of  parameters  and  sets  time  Intervals 
_ per  _plpt.« _  j 
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Fig  7B  Simulation  View  bum 
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For  simulation  of  manual  con- 

_ trol«_  either  an  air-to-ground  view 

_ from  on-board  the  payload  or  a 

_ ground-to-air  view  from  the  landing 

zone  can  be  selected.  A  variety  of  j 
three-dimensional  inertial  references, 

_ can_be_es tablished. . Limited  data  -  i 

_ perhaps,  only  elapsed  time  -  is  dis- I 

played  during  the  flight  with  full 
data  available  at  completion.  For 
simulation  of  automatic  guidance,  the 
appropriate  scheme  is  selected  and  j 
initialized  if  necessary.  Return  to  ! 
manual  control  will  occur  either  at  a 
selected  time  or  altitude.  A  running 
pictorial  display  Is  provided  to  j 

monitor  computer  performance  without  ; 
delaying  the  computation.  j 


Fig  s  Flight  Vata.  Unit 


Geometric  and  Aerodynamic 
Assumptions 


The  basic  components  of  the  pro¬ 
gram  r.re  the  assumed  aerodynamic  and  geometric  characteristics  and  the  | 
kinematic  relations.  The  canopy/payload  system  is  modeled  as  two  points  ! 
connected  by  an  in  extensible  line  subject  to  tensile  forces  (if  compressive 
forces  occur,  the  paylosd  f reef alls  until  the  orientation  is  such  as  to  j 
allow  tensile  forces) .  The  payload  mass  contributes  to  system  moments  of  I 
inertia  about  the  longitudinal  and  lateral  axes  only  since  suspension  geo-  j 
metry  usually  decouples  the  motion  about  the  vertical  axis.  A  constant  j 
jdrag  area  is  presently  assigned  to  the  payload  although  studies  of  particu-J 
lar  payload  configurations  with  lift,  drag,  and  side  forces  are  anticipated 
'The  mass  of  the  canopy  is  modeled  as  a  simple  rectangle  inclosing  an  air 
mass.  For  apparent  mass,  up  to  one-half  the  displaced  air  mass  (volume  of  : 
air  displaced  by  the  canopy  and  its  captive  internal  air)  is  selected  as  a  , 
constant  added  mass  concentrated  at  the  centroid  in  calculations  of  moments 
of  inertia.  The  geometric  distribution  of  the  captive  included  air  mass  is 
considered  in  calculation  of  moments  of  Inertia.  This  mass  aad  the 
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"  apparent  mass  .change. with  altitude  .according  to.  ambient  density.. _ ’ 

Certain  unique  aesumptions  have  been  made  in  the  treatment  of  aero- 
dynamic  forces  and  moments,  Bnpirical  functions  were  developed  for  the 
|basic  wing  lift  and  drag  coefficients  of  Parafoil  canopies  measured  in  the 
wind  tunnel.  A  side  force  ^efficient  varying  with  the  sine  of  the  side 
— h  ■  slip  angle  was  also  extraqtgd  from  the  wind  tunnel  Early  attempts  to! 

Incorporate  these  coefficients  and  other  moment  coefficients  and  stability  ; 
derivatives  measured  in  the  wind  tunnel  under  static  conditions  failed  to 
give  realistic  results.  Extremely  low.  turn. rates  were  noted  under  some 
typical  quasi-steady  conditions  even  .though  the .data  supposedly  reflected 

the  conditions  of  control  deflection. _ Realistic  performance  -  as  compared  ( 

with  motion  pictures  of  jumpers  and  limited  available  test  data  -  was 
finally  achieved  by  treating  the  canopy  as  geometrically  distributed  seg-  ; 

ments  Including  spanwlse  curvature  with. local. lif  t. jmd  drag^. computed ^or _ , 

each  segment  based  on.  local  velocity. and  .angle  of  attack.  Local  lift  co-  i 
efficients  vary  from  the  total  wing  lift  coefficients  in  a  manner  approxi-  1 
mating  an  elliptical  distribution.  The  empirical  functions  for  total  wing 
lift  and  drag  coefficients  Include  effects  of  trailing  edge  deflections.  i 
To  study  turn  and  pitch  (or  "braking")  control  with  trailing  edge  deflec¬ 
tion,  the  deflection  is  applied  via  the  empirical  functions  at  selected 
segments.  It  is  important  to  note  that,  due  to  the  spanwlse  curvature,  re¬ 
sultant  aerodynamic  forces  on  deflected  segments  may  be  significantly  in-  ; 
dined  to  the  vertical-longitudinal  plane. 


A  study  using  a  separate  computer  program  was  made  of  the  quasi-static 
stability  effects  of  typical  velocity  and  angle  of  attack  distributions 
resulting  from  the  segmentation  of  the  wing.  It  was  found  that  for  typical 
wing  geometries  associated  with  payloads  from  200  to  2000  pounds,  the 
effects  of  velocity  distribution  were  highly  stabilizing. (6) 


Experimental  Verification  of  Simulation 

I 

To  verify  the  simulation  program  a  flight  test  was  performed  using  a 
Parafoil  canopy  of  200  sq  ft  (18.6  sq  a)  and  a  payload  of  380  lbs  (165  Kg)  j 
specially  designed  to  measure  those  key  parameters  needed  for  verification.! 
The  emphasis  of  the  test  program  was  on  dynamic  flight  conditions  associated 
with  turning  maneuvers.  The  payload  Included  remote  radio  control  hard-  j 

ware  manually  activated  from  the  ground.  ! 

i 

Flight  Data  Unit  I 

i  * 

f  i 

S  The  flight  data  unit  (FDD)  serving  as  the  payload  consists  of  an  inner 
{box  cushioned  within  an  outer  shell  to  which  a  specially  shaped  foam  nose  ! 
cone  and  stabilizing  tail  fins  are  attached  (Fig  8).  The  geometry  of  the 
-dim;/  f ~*nR  snd  0096  cone  was  determined  during  wind  tunnel  tests  to  achieve  good  j 

•  in-  r  i  1  1 
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static  and  dynamic  stability.'  For  accurate  airspeed  and  gyro  data  the  pay- 
load  must  track  accurately  with  the  canopy.'  The  payload  must  als<r"track  ~ 
well  when  carried  by  sling  below  a  helicopter. 

A  J-Tec  vor t ex- sensing  airspeed  indicator  is  mounted  on  the  upper  part 
of  the  nose  cone,  Wind .tunnel  tests  were  made  for  calibration  at  various 
angles  of  attack.  Twin  ataxic  ports  in  the  surface  of  the  vertical  fins 
provide  static  pressure~£or  the  Rosemount  1241  altimeter  with  rate  output. 
This  mounting  arrangement  is  not  affected  by  moderate  alldesllp.  The  brand 
names  of  these  units  are  mentioned  because  they  appear  to  be  the  only  units 
available  that  are  adequate  for  this  application.  This  was  determined  in 
several  independent  surveys.  The  J-Tec  airspeed  indicator  acoustically 
senses  the  rate  of  vortex  shedding' 'from  a  special  cross-piece  in  the 
channel  of  the  indicator.  It  is  highly  accurate  and  responsive  in  the  test 
range  of  10  to  18  m/s.  Pressure  transducers  would  be  Inaccurate  in  this 
range.  The  Rosemount  altimeter  is  designed  for  test  applications  with  good 
accuracy  both  in  altitude  and  electronically-derived  rate  measurement. 
Airspeed  and  descent  rate  are  the  two  main  parameters  of  gliding  flight. 
Also,  the  horizontal  component  of  airspeed  can  be  calculated  from  the  vec¬ 
tor  difference  between  total  airspeed  and  descent  rate  only  if  each  para¬ 
meter  is  measured  accurately. 


The  FDU  also  contains  three  rate  gyros  on  the  yaw,  pitch  and  roll 
axes.  An  integrating  circuit  is  attached  to  the  yaw  rate  gyro  to  provide 
approximate  heading  data.  The  gyro  data  is  useful  in  indicating  the 
nature  of  turn  response  which  can  vary  widely  depending  upon  the  mechanics 
of  turn  control. 

Control  deflection  is  measured  from  the  servo  potentiometers. 
Originally,  current  draw  of  the  servo  motors  was  measured  as  an  indication 
of  control  force.  This  has  been  found  unreliable  as  an  Indication  of  force 
because  of  electrical  transients.  Modifications  are  being  made  to  measure 
forces  from  transducers  linked  to  the  lines. 


A  14-channel  FM  tape  recorder  is  used  to  record  all  signals.  A  play¬ 
back  unit  re-transforms  the  FM  signal  to  voltage  analog  to  drive  a  strip- 
chart  recorder  for  final  copy.  The  tape  recorder  is  satisfactory  for 
recording  data  during  a  flight  including  extreme  oscillations.  However, 
it  is  not  satisfactory  for  recording  data  under  high-shock  conditions  such 
as  during  deployment  and  landing.  Apparently,  the  tape  is  jiggled  as  it 
posses  the  head  when  the  unit  is  jolted. 

The  FDU  is  controlled  by  a  model  aircraft  radio-control  system.  The 
‘  radio-controlled  servos  drive  potentiometers  which  control  the  large  motors 
J  drawing  in  control  lines.  Currently,  two  channels  are  used  for  turn  and 
'  .  glide  ratio  control.  Three  motors  and  lines  are  used  for  the_left  turn. 
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right  turn  and  glide  ratio.  The  Parafoil  is  rigged  for  deflection  of  the 
third  (or  "C")  flare  on  the  second  set  In  from  each  tip  or  turn.  This  pro¬ 
duces  primarily  a  roll  tilt  of  the  wing  with  a  slight  Incidence  shift.  The 
trailing  edge  is  separately  deflected  for  glide  control. 

Weight  of  the  payload  jfas  355  pounds  during  initial  testing.  Later 
modifications  added  40  poia&s.  This  weight  is  carried  well  by  the  200- 
square-foot  Parafoil  with  unreefed  deployment  at  low  speeds  (10-30  knots) . 
Landing  speed  is  somewhat  high  with  this  loading. 

Except  for  the  airspeed  indicator,  all^ equipment  is  housed  within  the 
Inner  box  of  the  FDU  (35  X  22  X  14  Inches).  The  inner  box  is  tied  to  the 
outer  shell  at  the  bottom  and  surrounded  by  several  inches  of  polyurethane 
foam.  The  outer  shell  Is  108  inches  overall.  Including  nose  and  tail  with 
a  cross-section  of  35  X  28  Inches.  Drag  area  (CqS)  measured  approximately 
in  the  wind  tunnel  Is  2.15  -  2.69  square  feet  (0.2  -  0.25  sq  m)  . 

Setting  System  Parameters 

Those  aspects  of  the  simulation  model  needing  verification  Included: 
(1)  the  lift  and  drag  functions,  (2)  the  model  resulting  trim  condition 
predicted  by  the  longitudinal  stability  model,  (3)  the  assumption  that 
linearized  moment  coefficient  and  deviations  could  be  replaced  by  local 
force  computations  at  spanvise  locations,  and  (4)  the  computational  tech¬ 
niques  employed.  Thus,  in  matching  flight  data  none  of  these  factors  were 
altered.  Input  parameters  altered  for  the  match,  as  listed  below,  were 
those  parameters  not  physically  measurable.  Apparent  mass  (sometimes 
called  "virtual1*  mass)  is  a  flow  field  phenomenon  -  not  a  physical  mass. 
System  length  depends  on  In-flight  mass  distribution  and  distorted  suspen¬ 
sion  and  canopy  geometry.  While  these  parameters  could  be  "ballparked", 
they  could  not  be  physically  measured.  Ideally,  to  examine  flight  dynamics 
It  Is  desirable  to  start  from  a  steady  condition.  For  gliding  parachutes 
this  is  impossible.  After  over  twenty  flights,  steady  conditions  have 
never  been  observed.  Repeatedly,  a  maneuver  must  be  entered  before  the 
efforts  of  a  previous  maneuver  have  dissipated  even  though  flight  appears 
steady  from  the  ground.  Wind  turbulence  adds  additional  disruption.  The 
two  flights  matched  below  are  the  best  illustrations  for  verification  of 
the  simulation. 

To  obtain  a  reasonable  match,  the  following  system  parameters  were 
varied  in  the  simulation: 

1.  Apparent  mass  coefficient  (one  corresponds  to  half  the  volumetric 
displacement)  was  initially  assused  to  be  1.0.  Values  of  0,  0.8,  1.0,  1.5, 
2.0,  3.0  were  tried.  Final  value  la  1.0! 

2.  Fayload  drag  area  (CpS)  was  Initially  estimated  at  0.4  sq  m. 

Values  of  0.2,  0.25,  0.3,  0.35,  0.4  were  tried.  Value  selected  Is  0.2. 


GOOD  RICK 


Wind  tunnel  tests  later  confirmed  the  value  to  be  0.20  to  0.25  sq  m. 

3.  System  length,  which  Is  the  distance  from  payload  center~bf  mass 
to  canopy  center  of  mass,  was  varied  between  5.5  m  and  6,5  m.  Final  value 
is  6.5  m,  which  corresponds  with  physical  measurement. 

4.  Trailing  edge  deflection  was  assumed  and  confirmed  to  be  zero, 
though  deflection  of  the  edge  near  the  tip  was  assisted  in  some  trails. 

This  method  of  turn  control*  creates  high  initial  yaw  rates  which  were  not 
observed  in  the  flight  data. 

5.  Roll-tilt  turn  deflection  was  found  to  be  5.0  deg  in  order  to 
obtain  a  match.  Values  of  4.0  to  6.0  were  tested.  This  parameter  cannot 
be  measured  since  it  is  the  tilt  angle  of  the  resultant  force  vector.  The 
flight  data  Indicated  a  moderate  left  turn  at  neutral  control  line  deflec¬ 
tion  which  was  matched  by  assiming  -1.5  deg  (left)  tilt  at  neutral.  Thus, 
in  the  final  simulation  presented  here,  tilt  varies  from  -1.5  to  5.0  deg 
for  1002  deflection. 

6.  Incidence  coupling  was  found  to  occur  with  tilt-turn  deflection. 
This  was  the  main  surprise  learned  in  the  matching  process.  The  incidence 
angle  is  the  angle  of  the  chord  below  the  canopy  X-axis  with  the  Z-axis 
passing  from  center  of  mass  to  payload  center  of  mass.  In  the  stability 
analysis  of  Reference  3,  the  optimum  value  was  found  to  be  4  deg  which 
appears  close  to  the  constructed  value.  Film  of  the  flight  test  showed 
substantial  incidence  change  of  the  section  to  which  the  control  line 
attached.  Also,  it  was  impossible  to  match  the  characteristic  Initial  de¬ 
crease  in  descent  rate  as  turn  deflection  was  applied  if  incidence  angle 
was  held  constant  at  4  deg.  But  when  incidence  was  allowed  to  decrease 

to  2  deg  linearly  as  turn  deflection  went  to  1002,  the  proper  decrease  in 
descent  rate  was  obtained,  and  the  proper  sensitivity  toward  spin  was 
observed.  Though  no  spin  occurred  in  the  first  flight,  considerable  sensi¬ 
tivity  to  duration  of  the  applied  stroke  was  observed.  Thus,  to  obtain  the 
data  match.  Incidence  was  assumed  set  at  4  deg  (nose  down)  varying  to  2  deg 
(nose  down)  with  1002  turn  deflection. 

Comparison  of  Data 


The  descent  rate  data  matches  well 
in  parts  with  certain  deviations  (Fig 
9).  Descent  rate  is  the  key  parameter 
for  comparison  since  all  force,  moment 
and  orientation  variations  have  a  dir¬ 
ect  effect  on  it.  Mote  the  large  de¬ 
viation  from  equilibrium  (V2  “  -4.7 
m/s) .  This  indicates  that  the  turns 
are  highly  unsteady  due  to  excessive  de¬ 
flection.  Turns  of  moderate  rate  would 
show  excursions  to  -6  or  -7  m/s  after 
deflection  with  little  sensitivity  to 
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duration  of  the  deflection.  In  the  simulation  just  previous  .to  the  one  pre¬ 
sented,  the  second  turn  deflection  was  a  few  percent  greater  initially, 
producing  a  peak  descent  rate  of  -16  m/s.  If  the  first  turn  were  held  one 
to  two  seconds  longer,  a  high  rate  spin  would  occur  with  descent  rate  peak¬ 
ing  at  -22  m/s  and  a  turn  rate  of  ISO  deg/sec.  A  spin  almost  that  bad 
occurred  on  a  later  flights 


Fig  10  Total  and  Horizontal  KUuptzd 


Fig  11  You)  Rate 


Fig  10  shows  the  total  airspeed 
measured  and  the  horizontal  airspeed 
computed  as  the  vector  difference  of 
total  with  descent  rate.  It  was 
difficult  to  match  initial  conditions 
since  the  system  was  just  recovering 
from  a  previous  right  turn.  It  is 
encouraging  that  horizontal  airspeed 
can  be  accurately  estimated  from 
measurement  of  total  airspeed  and 
descent  rate  because  this  allows 
evaluation  of  glide  performance  in¬ 
dependent  of  wind  and  also  allows 
estimation  of  wind  speed  when  com¬ 
pared  to  measured  speed  by  a  guid¬ 
ance  system  (not  done  for  these 
flights) . 

Fig  11  shows  the  axial  yaw  rate 
data.  Fig  12  shows  yaw,  pitch  and 
roll  rates  for  the  first  turn.  The 
gyro  data  includes  significant 
"noise”  caused  by  payload  oscilla¬ 
tions.  The  payload  oscillations  are 
of  high  rate  (10  deg/sec)  but  of 
short  period  so  that  system  motion 
can  be  identified.  It  was  not  known 
whether  gyro  data  would  be  meaning¬ 
ful  or  not.  For  turning  motions 
encountered  in  this  flight  there  is 
reasonable  correlation  at  the  peaks 
on  all  axes. 


Spin  Performance 


An  unintended  spin  produced  data  which  was  also  matched  closely  by 
simulation.  During  deployment,  a  suspension  line  wrapped  over  one  cell  at 
the  tip,  distorting  the  outer  25Z  of  the  canopy.  Measured  descent  rate 
averaged  9.7  m/s  with  oscillations  of  between  0.5  and  3  m/s  with  a  period 
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of  1,6  seconds,  Film  indicated 
that  the  trailing  edge  was 
partially  deflected  as  a  result 
of  twisted  lines.  Otherwise, 
most  of  the  canopy  was  Inflated 
to  a  proper  shape.  Simulations 
were  made  using  zero,  half,  and 
full  flap.  Descent  rates  aver¬ 
aged  22,  10,  and  8  m/s,  respec¬ 
tively.  At  half  flap  the  simula¬ 
tion  also  showed  (Fig  13)  an 
oscillation  of  2  m/s,  about  the 
average  of  10  m/s  damping  to  0.5 
m/s  in  11  seconds.  Period  of 
oscillation  was  1.6  seconds.  The 
Euler  (or  inertial)  yaw  rate  dur¬ 
ing  the  flight  was  estimated  from 
the  video  tape  at  197  deg/sec 
(averaged  over  ten  revolutions) . 
Simulation  at  half  flap  indicated 
205  deg/sec.  Other  data  normally 
recorded  was  off-scale  during 
this  flight.  The  close  match  of 
rate  of  descent  and  yaw  rate  adds 
additional  validity  to  the  simu¬ 
lation  model.  The  fact  that  this 
occurred  in  spite  of  the  canopy 
distortion  may  be  explained  by 
the  assumption  that  the  distorted 
portion  would  be  stalled  anyway 
in  a  spin,  producing  only  drag. 
Spin  dynamics  appear  to  be  driven 
mainly  by  the  mass  distribution. 
The  flight  data  did  not  exhibit 
the  same  damping  characteristics; 
however,  the  system  may  have 
. .  been  subject  to  various  excita¬ 
tions  as  the  payload  pitched, 

wind  varied,  or  other  changes  occurred.  The  simulation  was  controlled  by 
holding  full  left  turn  and  half  flap  while  starting  at  a  nominal,  sero- 
flap  glide  condition.  The  spin  was  established  within  3  seconds.  Curious¬ 
ly,  the  ground  track  showed  tight  left-turning  circles  whose  centers  pre- 
cess  circularly  to  the  right. 

Comparison  has  been  made  in  a  few  other  cases  with  results  similar  to 
those  presented  here.  Therefore,  the  experimental  results  appear  to  vali- 


fig  13  Descent  Rate  Voting  Spin 
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date  the  simulation  program  as  applied  to  a  200-square  foot  Parafoil  of 
aspect  ratio  of  2.0  with  a  350-pound  payload.  This  conclusion  is  warranted 
because: 

1.  Due  to  the  discrete  character  of  control  inputs,  response  is 
driven  mainly  by  aerodynamic  and  Internal  properties  assumed  for  the  model. 

2.  The  combination  oL*ulti-axial  responses  predicted  by  the  model 

Is  in  close  agreement  with  Experimental  results. _ 


$ 

0 

, 

'  x 


As  of  this  writing,  about  10  good  data  flights  have  been  made.  The 
simulation  model  has  matched  salient  features  of  most  flights.  There  is 
some  difficulty  when  wind  changes  occur.  The  data  shows  a  response  to  a 
disturbance  rather  than  a  control_  input  without  Indicating  sufficient  data 
for  simulation.  This  has  made  it  difficult  to  match  flare  dynamics  since 
/  the  few  flare  landings  made  were  accompanied  by  some  wind  turbulence.  The 
simulation  has  been  used  to  identify  some  anomalies  such  as  incidence 
vshifts  due  to  line  stretch  and_ trim  change  caused  by  a  pilot  chute. 


\  Conclusion 

\  \ 

The  extent  of  agreement  between  simulation  and  flight  data  indicates 
that  the  primary  factors  included  in  the  longitudinal  stability  analysis 
and  in  the  6DOF  simulation  are  correct.  Although  the  stability  analysis 
predicts  only  steady-state  behavior,  it  forms  the  basis  required  for 
analysis  of  dynamic  behavior  in  the  body-fixed  XZ  plane.  The  agreement 
seen  in  descent  rate  ^(FS$-9>  would  not  be  possible  with  an  invalid  sta¬ 
bility  model.  Of  course,  in  turning  flight  during  the  response  Immediately 
following  the  deflection,  other  factors  such  as  the  assumed  spanwise  distri¬ 
bution  of  lift  and  drag  become  predominant.  The  agreement  in  yaw  rate 
(Pi^lll  best  illustrates  correctness  of  this  aspect  of  the  6D0F  model. 

The  mechanics  driving  motion  during  a  spin  are  quite  difficult  to  under¬ 
stand;  however,  the  agreement  shown  in  the  descent  rate’^tF^— 1-3)  indicates 
that  the  mass  ratios  assisned  are  accurate  and  further  justifies  the  assump¬ 
tion  of  spanwise  distribution  of  lift  and  drag. 


In  further  development  activities  on  gliding  parachute  systems,  the 
6D0F  model  will  serve  to  guide  exploratory  work  and  will  be  updated  for 
more  accurate  application  to  different  canopies  and  to  larger  systems., 

/ 

References 

\ 

1.  Goodrich,  T.,  "Wind  Effects  on  Gliding  Parachute  Systems  with  Non- 
pro  port  ional  Automatic  Homing  Control",  TR70-28-AD,  USA  Natick  R&D  Labora- 

,  tor lea,  1969. 

2.  Goodrich,  T.,  "Estimation  of  Wind  Effect  on  Gliding  Parachute 
.Cargo  Systems  Using  Computer  Simulation",  Paper  #70-1193,  American  Insti¬ 
tute  for  Aeronautics  and  Astronautics  (AIAA),  1970. 


I 

I 


i 


f  j 


ra  r  i  <-•: 


i 


GOODRICK 


3.  Good  rick,  T,_,_Pearson,  A.,  Murphy,  A.,"  Ana  lysis  ofVarious  Auto¬ 
matic  Honing  Techniques  for  Gliding  Airdrop  Systems  with  Comparative 
Performance  in  Adverse  Winds",  Paper  #73-462,  American  Institute  for 
Aeronautics  and  Astronautics  (AIAA),  1973. 

4.  Ware,  George  M.  and  Hassell,  Jones  L.  Jr.,  "Wind  Tunnel  Investl~ 

gatlon  of  Ram-Air-Inflated  All -Flexible  Wings  of  Aspect  Ratios  1.0  to  3.0", 
NASA  1M  SX-1923,  Langly  Res&rch  Center,  1969,  ..  . . 

5.  Goodrich,  T.,  "Theoretical  Study  of  the  Longitudinal  Stability  of 
High-Performance  Gliding  Airdrop  Systems",  Paper  #75-1394,  American  Insti¬ 
tute  for  Aeronautics  and  Astronautics  (AIAA),  1973. 

6.  Goodrich,  T.,  "Simulation .Studies  of  Jthe  Flight  Dynamics  of  Glid-  ' 
ing  Parachute  Systems",  Paper  #79-0417,  American  Institute  for  Aeronautics 
and  Astronautics  (AIAA),  1979. 


